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Abstract 

This paper deals with feedback control of fractional-order Chua's system. 
The fractional-order Chua's system with total order less than three which 
exhibit chaos as well as other nonlinear behavior and theory for control of 
chaotic systems using sampled data are presented. Numerical experimental 
example is shown to verify the theoretical results. 



1 Introduction 

It is well known that chaos cannot occur in continuous systems of total order 
less than three. This assertion is based on the usual concepts of order, such as 
the number of states in a system or the total number of separate differentiations 
or integrations in the system. The model of system can be rearranged to three 
single differential equations, where one of the equations contains the non-integer 
(fractional) order derivative. The total order of system is changed from 3 to 2-|-g, 
where < g < 1. To put this fact into context, we can consider the fractional- 
order dynamical model of the system. Hartley et al. consider the fractional- 
order Chua's system and demonstrated that chaos is possible for systems where 
the order is less than three. In their work, the limits on the mathematical order of 
the system to have a chaotic response, as measured from the bifurcation diagrams, 
are approximately from 2.5 to 3.8. In work chaos was discovered in fractional- 
order two-cell cellular neural networks and also in work ^ chaos was exhibited 
in a system with total order less than three. 

The control of chaos has been studied and observed in experiments (e.g. works 
0, ||12|). Especially, the control of well-known Chua's system by 



sampled data has been studied [|l^. The main motivation for the control of 
chaos via sampled data is well-developed digital control techniques. 

In this brief studies are presented practical results from sampled-data feed- 
back control of fractional-order chaotic dynamical system, modeled by the state 
equation x = f(x), where x G 3ft" is state variable, f : 3ft" — »• 3ft" is nonlinear 
function and f(0) = 0. 

The approach used in this paper is concentrating on the feed-back control of 
the chaotic fractional-order Chua's system, where total order of the system is 2.9. 



2 Fractional calculus 



2.1 Definitions of Fractional Derivatives 

The idea of fractional calculus has been known since the development of the 
regular calculus, with the first reference probably being associated with Leibniz 
and L'Hospital in 1695. 

Fractional calculus is a generalization of integration and differentiation to 
non- integer order fundamental operator a-D", where a and t are the limits of the 
operation. The continuous integro-differential operator is defined as 

r 3f?(«) > 0, 

aDt = { 1 3?(«) = 0, 

[ /,*(rfr)-" 3fJ(a) < 0. 

The two definitions used for the general fractional differintegral are the Griinwald- 
Letnikov (GL) definition and the Riemann-Liouville (RL) definition 0, The 
GL is given here 

^D'^m = \imh--Y.i-^y( Mit-Jh), (1) 
where [x] means the integer part of x. The RL definition is given as 

for in — \ < a <n) and where Ff.) is the Gamma function. 



2.2 Numerical Methods for Calculation of Fractional 
Derivatives 

For numerical calculation of fractional-order derivation we can use the relation 
derived from the Griinwald-Letnikov definition (p. This approach is based 
on the fact that a wide class of functions, two definitions - GL (|1]) and RL (Q) 
- are equivalent. The relation for the explicit numerical approximation of a-th 



derivative at the points fcT, (fc = 1, 2, . . .) has the following form 0, |T0|, |Tl[| : 

ik-L/T)D^,Tfit) ^ T-"x:(-i)^f"V^-^' 

j=0 \J ) 

where L is the "memory length", T is the step size of the calculation (sample 
period) and (—!)■' (^"^ are binomial coefficients Cj"'*, (j = 0, 1, . . .). For its calcu- 
lation we can use the following expression: 



1 + a 



= 1, of =1-^1 cr_v (4) 




2.3 Some Properties of Fractional Derivatives 



Two general properties of fractional derivative we will be used. The first is com- 
position of fractional with integer-order derivative and the second is the property 
of linearity. 



The fractional-order derivative commutes with integer-order derivation [11 



^iaDUit)) = aDl (^^) = (5) 

under the condition t = a we have f^^\a) = 0, (A; = 0, 1, 2, . . . , n - 1). The 
relationship says the operators and commute. 

Similar to integer-order differentiation, fractional differentiation is a linear 



operation 1 1 



„Df (A/(t) = XaDU{t)+iiMg(t). (6) 



3 Fractional- Order Chua's System 

Classical Chua's oscillator, which is shown in Fig. 1, is given by 

It [G{v, - V2) + , (7) 

-'2 

\-[-V2{t)-RLz], 

where G = 1/R and /(fi) is the piecewise linear v — i characteristic of nonlinear 
Chua's diode. 

Given the techniques of fractional calculus, there are still a number of ways in 
which the order of system could be amended. One approach would be to change 
the order of any or all of tree constitutive equations (0) so that the total order 
gave the desired value. 

In our case, in the equation one, we replace the first differentiation by frac- 
tional differentiation of order g, g G 3?. The final dimensionless equations of the 
system for = are (xi = Vi,X2 = f2, = i): 

a qD]^" {x2it) - xi{t) - fixi)) , 

Xi{t)-X2it)+X3{t), (8) 
-Px2{t), 
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Figure 1: Chua's circuit. 

where 

/(xi) = bxi + -{a - b){\xi + 1| - \xi - 1|) 
and a = C2/C1, f3 = C2R^/L. 

4 Feedback Control of Chaos 



The structure of control system with sampled data [T^ is shown in Fig. 2. 



The state variables of the chaotic system are measured and the result is used to 
construct the output signal y(t) = i5x(t), where D is a constant matrix. The 
output y(t) is then sampled by sampling block to obtain y(fc) = Dy.{k) at the 
discrete moments /cT, where k = 0, 1, 2, . . ., and T is the sample period. Then 
D^[k) is used by the controller to calculate the control signal u(fc), which is fed 
back into chaotic system. 



The controlled chaotic system is defined by relations |13] 

= f(x(t)) + 5u(A;), te[kT,{k + l)T) 

u(fc + l) = Cu(A;) + Dx(A;), A; = 0,1,2,... (9) 

where ue3?'", 5 G 3fJ''x3fJ'", C e a^^xa?"", D e afJ'^xgfJ" and t G x(fc) is the 
sampled value of x(t) at t = kT . Observe that since f(0) = is an equilibrium 
point of the system (^. 

The controlled fractional-order Chua's system is defined by 

dxi{t) ^ 



dt 
dx2{t) 

dt 
dxsjt) 

dt 



a oA"^ ix2{t) - xi{t) - f{xi)) + ui{t), 

Xiit)-X2it)+X3{t)+U2it), (10) 
-pX2{t)+Us{t). 
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Figure 2: Structure of control system. 

5 Illustrative Example 

For numerical simulations the following parameters of the fractional Chua's sys- 
tem (H) were chosen: 

a = 10, p = q = 0.9, a = -1.27, b = -0.68, 
and the following parameters (experimentally found) of controller: 
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Using the above parameters (IT^) the digital controller in state space form is 
defined as 

ui{k + l) = 0.8ui{k) -3.3xi{k), (12) 

for A; = 0, 1, 2, . . . The initial conditions for Chua's circuit were ((a;i(0), ^2(0), ^3(0)) 
= (0.2, —0.1, —0.01) and the initial condition for the controller (p!2D was {{ui{0) 
= (0)). The sampling period (frequency) was T = 100 Hz. 

For the computation of the fractional-order derivative in equations (^^), the 
relations (^, (^ and properties (^, (^ were used. The length of memory was 
L = 10 (1000 coefficients for T = 100 Hz). 

Fig. 3 shows the attractor of Chua's circuit (^ without control. Similar 
behaviour was shown in work |^ , where piecewise linear nonlinearity was replaced 
by cubic nonlinearity which yields very similar behaviour. 



Figure 3: Strange attractor of fractional- order Chua's system without control. 




Figure 4: Controlled state variables Xi{t),X2{t),x^{t). 



Figure 5: Control signal Ui{t). 



In Fig. 4 is shown the controlled trajectory of state variables of the fractional- 
order Chua's system (p!OD, which tends to origin asymptotically. 
In Fig. 5 is shown control signal from the digital controller ([T^). 



6 Conclusion 

We have considered an example of control of chaotic fractional-order Chua's cir- 
cuit, which exhibits chaotic behaviour with total order less than three. As has 
been demonstrated, the idea of fractional calculus requires one to reconsider 
dynamic system concepts that are often taken for granted. So by decreasing 
the order of a system from 3 to 2 + g in this way, we also move from a three- 
dimensional system to one of infinite dimension. This system can be controlled 
by sampled data. The sampled data of output are sufficient for constructing the 
control signals in the digital controller. Digital controllers had been widely used 
in industry. 



The conclusion of this work confirms the conclusions of the works []TT 
that there is a need to refine the notion of the order of a system which can not 
be considered only by the total number of differentiation. For fractional-order 
differential equations the number of terms is more important than the order of 
differentiation. 

The results presented in this contribution give basis for controlling chaotic 
fractional-order systems. An alternative approximation of fractional-order deriva- 
tive, stability investigation, and also other chaotic fractional-order system will be 
used in further work. 
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